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Glide and climb of quantum dislocations under finite external stress, variation of chemical po- 
tential and bias (geometrical slanting) in Peierls potential are studied by Monte Carlo simulations 
of the effective string model. We treat on unified ground quantum effects at finite temperatures 
T. Climb at low T is assisted by superflow along dislocation core - superclimb. Above some crit- 
ical stress avalanche- type creation of kinks is found. It is characterized by hysteretic behavior at 
low T. At finite biases gliding dislocation remains rough even at lowest T - the behavior opposite 
to non-slanted dislocations. In contrast to glide, superclimb is characterized by quantum smooth 
state at low temperatures even for finite bias. In some intermediate T-range giant values of the 
compressibility as well as non-Luttinger type behavior of the core superfluid are observed. 
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I. INTRODUCTION 

Superfluid and mechanical properties of solid 4 He continue to be fascinating and, in many respects, 
largely enigmatic subjects for study. This is especially so for the proposed supersolid phase — and 
the observed anomaly interpreted as non-classical rotational inertia (NCR!) in torsional oscillator (TO) 
experiments — ^. First-principle simulations — as well as a strong dependence of the NCRI on preparation 
and annealing protocols (see in Ref£) have definitely excluded supersolidity of ideal hep solid 4 He. Thus, 
the superfluid response of solid 4 Hc can only be due to superfluidity along extended structural defects 
forming a percolating network. Dislocations are one of the possible candidates for such network - some 
of them do have superfluid cores as determined in ab initio simulations^—. 

Dislocations control mechanical properties of solids. One of the remarkable effects - temperature 
induced shear modulus softening^ has recently been observed in solid 4 He as well^. There is the whole 
plethora of very intriguing and largely unexplained features exhibited by the modulus at extremely 
small external stresses^. Very recently the opposite behavior - low-T softening has been reported^ 1 - for 
impurities free crystals. 

Superclimb of edge dislocations^ is macroscopic quantum phenomenon where high-T diffusion along 
dislocation core is replaced by superfluid transport. Convincing evidence for such transport has been 
presented in RefsJ^. There are three crucial aspects of the observation of mass transport through solid 
4 Hei£: i) Suppression above temperature To « 0.6K which is significantly lower than that for the bulk 
A point; ii) Uniform transfer of time dependent pressure through the whole sample during the superflow 
events ; iii) Macroscopically large mass accumulation in the bulk during the superflow events. The feature 
i) clearly excludes any macroscopic superfluid channels scenario. The item ii) is inconsistent with classical 
plasticity driven by gradients of stress. The effect iii) finds its natural explanation in terms of superfluid 
mass transport along edge dislocation cores which adds (or removes) material to (from) half-plane forming 
edge dislocations^. 

In this paper, we analyze both aspects of dislocation thermodynamics - their glide and superclimb. 
Specifically, we consider stress induced roughening of gliding dislocation in Peierls potential. Finite 
biases due to pinning of dislocation ends in different valleys of Peierls potential are also analyzed. We 
also consider temperature dependence of so called isochoric compressibility x(T) and specifically address 
the prediction of its suppression at low T made in Ref4. We propose utilizing two types of responses in 
experimental studies full and differential. 



II. STRING MODEL OF DISLOCATION 
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Wc model dislocation as a quantum string of length L x strongly pinned at its both ends and subjected 
to Pcicrls potential Up, external force F and slanting (bias). Such string can perform two types of motion 
- glide and superclimb. The model for glide has been discussed in detail in Ref^. The modification 
of this model to describe superclimb has been introduced in Ref.— . Here we present unified description 
suitable for both cases, beginning with the superclimbing dislocation. After trivial simplification, the 
resulting action can be utilized for glide as well. In this paper we do not consider more involved model - 
dislocation performing glide and superclimb. 

Thcrmodynamical properties of dislocation arc described by its partition function Z = 
J Dy(x, t)Dp(x, t)D(j)(x, t) exp(-H) calculated in imaginary time < t < ft = 1/T (units H = 1, Kp = 1), 
where the functional integral is evaluated with respect to the dislocation displacement y(x,t), density of 
the core p(x, t) and superfluid phase 4>(x, t) conjugate to p, with x being coordinate along the dislocation. 
The action 

H = Hsf + Hdis (1) 

H SF = f dt^iP + no^t^+^^ + ^-ip-y) 2 } (2) 
Jo x 2 2p 

H dis = /' /3 ^[-^-((V t2 ;) 2 + y d 2 (V x2 /) 2 )- U pCos(27ry) 

JO x Z1S -d 

+ ^y 2 -Fy(x,t)] (3) 
y(0,t) = 0; y(L x ,t) = y max = 0,1,2,..., y(x,t + 0) = y(x,t). (4) 

is in continuous time t and discrete space x = 0, 1, 2, ...L x . Here Hsf describes core superfluidity^, with 
no,po standing for average filling factor and bare superfluid stiffness, respectively; units are such that 
speed of sound in superfluid is unity. In this work we are considering the case of integer no- Non- integer 
no lead to phase separation effects at low T and require additional analysis which will be presented 
elsewhere. 

The derivative V x (f> is understood as a finite difference defined modulo 2n: that is, V x <j)(x,i) = 
cos((f>(x + 1, t) — 4>(x, t)). In what follows, cos will be treated in the Villain approximation cos(9(x, t)) — > 
9(x,t) + 2irm x (x,t), with the ensemble summation over the integer variable m x {x,t) — 0, ±1, ±2, ... 
performed. 

We discretize time into N t slices and replace Vt0(x,t) —> [4>(x,t + 1) — 4>{x,t) + 2nmt(x, £)]t -1 , with 
m t (x,t) = 0, ±1,±2, and r standing for the time slice chosen as t = l/(TN t ). The last term 
in Eq.Q describes the superclimb effect - building the edge of the dislocation so that its no energy- 
cost displacement y — > y ± 1,±2, ... is possible by delivering matter p — > p ± 1, ±2, respectively, 
through supcrflow along the core. All the variables are periodic along the imaginary time /3. We consider 
dislocation attached to large superfluid reservoirs at its both ends. This implies periodic in space boundary 
conditions for supcrcurrents. 

The term ([3]) describes string motion, with l/Kd,Vd standing for dislocation effective mass density 
and speed of sound (in units of speed of sound in superfluid), respectively, and u p denotes strength of 
Peierls potential, with Burgers vector b taken as unity. The term ~ uo accounts for effects of dislocation 
pinning by 3 He impurities (not discussed here) and various structural defects. The last two terms in 
Eq.Q account for the effect of external stress a so that the linear force density F w ba (where we ignore 
spatial indices). The last line, Eq.(Q|, describes strong pinning at x = and x = L x so that, in general, 
the string is biased by y m ax ^ to have finite number of kinks for gliding dislocation (and - jogs for 
superclimbing) even at T = 0. 

We treat the Peierls term in the Villain approximation as well, that is, — cos (2iry(x, t)) — >~ (y(x, t) + 
P(x,t)) 2 with the ensemble summation over the discrete variable P(x,t) = 0, ±1,±2, .... This allows 
exact integrating out the variables p, <p, y analytically and expressing the partition function solely in 
terms of the integer variables m x (x, t), mt{x, t), P(x, t). Integration of the phase generates the constraint 
Vtm t + V x m x — 0, where now V*, are understood as simple discrete differences with respect to the 



integer variables m x ,mt denned on the bonds of the space-time lattice. Such constraint constitutes so" 
called J-current modelii which was treated within the Worm Algorithm (WA}». The variable P is 
unconstrained. Monte Carlo simulations have been performed for fixed number of time slices N t (we 
used N t = 40 in this work) with updates typical for WA modified in the presence of independent integer 
variable P defined on each site of the space-time lattice. Discrete time simulations carry systematic 
error decreasing with increasing Nt (for fixed T). However, we found that thermal profiles of measured 
quantities are self-similar for different Nt > 20 and, therefore, such error can be absorbed into the 
definition of unit of temperature. The value Nt — 40 is a reasonable choice between acceptable systematic 
error and realistic simulation times. 



A. Long-range interaction 

The parameter Kd in Eq.Q is not actually a constant. As discussed in RefJ^, it contains a contribution 
due to long-range forces between kinks. At long distance r between two kinks the interaction potential is 
~ 1/r. Accordingly, 1/Kd has a log-divergent facto r 13 ' 16 . Ignoring retardation effects (which do not lead 
to any significant differences), in Fourier 

K^iq) = ±C{q), C(q) = l + U c ln(l + j^—\, (5) 

where K ~ 1, JJ C ~ 1 are parameters; b ~ 1 stands for the short scale cutoff; q is spatial Fourier 
vector along the dislocation. As discussed in Ref.— , in the absence of the long-range forces, Uc = 0, the 
parameter K determines if the system is quantum rough (for K > K c with K c = 2/n) or quantum smooth 
(K < K c ). At K = K c it undergoes Berezinskii-Kosterlitz-Thouless (BKT) transition at T = driven 
by quantum fluctuations (for long strings). In the quantum rough state, Peierls potential is irrelevant 
and dislocation behaves as free Granato-Lucke string^. One can estimate (in the standard units up to 
a numerical factor ~ 1) K = Ha/(moVdb 2 ), where uiq is atomic mass, a stands for interatomic distance 
along the dislocation core; b - Burgers vector; Vd - speed of sound in standard units. Evaluation of K 
for solid 4 He, shows that K ~ 1, that is, the situation is marginal (as compared to crystals of heavier 
elements or those with significantly higher Vd where K « 1). However, as proven in Refi^ 3 -, arbitrary 
small long-range interaction (Jp^O always keeps unbiased dislocation in the smooth state at T = 0. In 
contrast to glide, superclimb is not affected qualitatively by such long-range interaction^, and dislocation 
becomes quantum smooth at low T for any value of K even for Uc = 0. 



III. GLIDE 



In this section we will discuss gliding dislocation at finite external stress and bias. The above model 
can be reduced to describe glide as considered in Refj^ 3 - by setting p s — > oo so that Hsf can be simply 
set to zero. 

We introduce two responses, one is differential 

R2 = J_ ^ d(y(x,t)) _ 1 d 2 In Z(F) 
pL 

and another one (we call it full) 



f ^ d{y(x,t)) 1 In Z(F) 

Jo ^ dF ~ (3L X dF* W 



1 ft (y(x,t)) 1 d\nZ(F) 



f3L x J ^ F pL x FdF 

where the averages are taken over the full action ([TJ with finite F and bias y m ax- Since at finite biases 
there is already finite (y) ~ y m ax, we modify Eq.([7]) by requiring that R\ includes average of two cases 
±Umax- This requirement reflects a generic situation that actual samples contain on average equal amount 



of oppositely slanted dislocations. In what follows we will see that, while at small stresses these two 
responses coincide, at finite stresses they behave quite differently, and, therefore, important information 
can be inferred about state of dislocations. Physically, two responses can be measured at slightly different 
frequencies so that total applied stress a = <tq + a' consists of finite part <7q at one frequency and small 
addition a' - at another. Then, i?2 can be viewed as linear response on a' at finite a§. 

i?i : 2 can be related to, respectively, full G\ and differential G2 shear moduli following the procedure 
described in detail in Refill. If the dislocation network is viewed as a set of blocks of free segments of 
sizes L x ,L y ,L z along the corresponding orthogonal axes, a displacement y (along y-axis) under a force 
F = <7 zy b where a zy stands for the zy-component of the stress tensor, results in a strain of the block 
u zy ps by/(L y L z ) with y oc a zy . Specifically: 

u zy = ^ [ dtJ2(vM, N d = -±-, (8) 

P^x Jo ~ ^y^z 

where Nd can be called dislocation density with respect to chosen coordinates. We introduce dimensionless 
dislocation density rid = b 2 Nd in units of b. Clearly, our description is based on the assumption rid << 1- 
It should also be noted that, in reality, dislocation density is not a scalar. We, however, will be working 
within the simplifying assumption L y ps L z ps L x , so that rid ~ b 2 /L 2 can be treated as a scalar for 

all practical purposes. Eq.lJHJ describes macroscopic plastic deformation. Adding elastic strain u> zy = 
G~ 1 a zy , where G e denotes bare elastic shear modulus, and relating full strain to stress, one finds 

1111 

77- = yr + n d Ri, yr = Tr + n dR2, (9) 

(_rl Lr e Lr2 <Je 

At high T Peicrls potential becomes irrelevant, and the responses become equal to each other and 
to that of free string R ps (1/L x ) Y^ x 0.5KV^~ x(L x — x) ps KV d ~ 2 L 2 /12, where we ignored log-factor 
([5|). Thus, given G e KV^ ps 1 and that ridL 2 ps 1 for isotropic dislocation forest, one concludes that 
high-T values Gioo = G200 = G^ ss G e /(1 + 7) where 7 ps 0.1 is a geometrical factor independent of 
dislocation density. Should, however, the dislocation forest has some asymmetry, say, L 2 >> L y L z , the 
factor 7 ps 0.1L 2 / L y Lz can, in principle, take arbitrary large value. It is convenient to introduce i?i ; 2 
relative to Rq for given L x . So, in what follows we will be discussing namely such values R12 = Ri,2/Ro, 
approaching unity at high T. Then, Eqs.Q can be rewritten as 

Gi l2 - Ge ~ (10) 

1+7^1,2 



As proven in RefJ^, at T = gliding non-slanted dislocation must be pinned by Peierls potential. 
This implies R\ t 2 nd << 1 so that G\^ — G e . As T increases thermal roughening takes place due to 
creation of kink-antikink pairs. This leads to increase of J?i,2 so that eventually R1.2 1 and the moduli 
G\.2 attain their high-T value Goo. RcfJ^has focused on describing this effect in the limit of no bias 
{Umax = 0) and vanishing external stress (F — > 0). 

Here wc address the question how finite stress and bias affect the roughening process. In general, 
roughening could be expected to be achieved by increasing temperature, bias and applied stress. All 
three factors tend to introduce kinks which may be expected to behave as free quantum particles which 
liberate dislocation from Peierls potential. However, as we will discuss below, these three factors lead to 
different scenarios. 



A. Finite stress 



Thermal roughening at small stresses is a relatively smooth crossover—. It is demonstrated by lower 
curve in FigfTJ In contrast, inducing roughening by increasing stress is characterized by non-monotonous, 
jump-like features. Fig J5] demonstrates such behavior for Ri and R2 at constant temperature T = 0.02 (in 
dimensionless units). It is important that the threshold stress (the parameter F), inducing such features, 
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FIG. 1: (Color online) R2 linear (F = 0.001) response of non-slanted (lower curve) and slanted (upper curve) 
dislocations. Parameters are L x — 80, Uc = 2, K = 1, Va = 1; bias is determined as linear density of extra kinks. 
Inset shows corresponding G2-modulus determined from Eq. (|10[) with 7 = 0.15. Slanted dislocations induce low- 
T shear modulus softening shown by the lower curve in the inset. Green gashed line indicates response of free 
Granato-Liicke string^. 
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FIG. 2: (Color online) Full Ri and differential R2 responses of gliding dislocations of different lengths on applied 
external stress a (F — ba) at temperature T = 0.02. Green dashed lines indicate responses of free strings (in the 
absence of Peierls potential). Error bars are shown for every point. Lines are guides to eyes. Note that the 
threshold F c for F (sharp increase of Ri and peak in R2) is about two times smaller for L x — 80 than that for 
L x = 40. Next peaks positions are approximately periodic with the period ~ F c . 
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FIG. 3: (Color online) Energy histogram for gliding dislocation undergoing low-T stress induced roughening. 
Bimodality of P(E) is a clear signature of jump-like transition and hysteretic behavior. 

depends essentially on the dislocation length. The non-linearity is associated with work FL X done by 
stress to create a single kink-antikink pair carrying energy AiS. Further increase of F creates more kinks 
in a quasi-periodic fashion with the overall tendency to approach free string response (green dashed line) . 
This behavior is clearly seen for L x = 80 case in Fig|2] 

We have found that force induced roughening proceeds similarly to I-order-transition at low T. Fig[3] 
demonstrates bimodal distribution of energy for dislocation undergoing stress induced roughening at 
T = 0.01. Clearly, such bimodality implies hysteresis in mechanical properties. The bimodality vanishes 
above certain temperature. For chosen parameters Uc = 2,u p = 0.1, K = 1,L X = 80 this occurs 
above T ss 0.015. More systematic studies how the onset temperature of hysteretic behavior depends on 
dislocation parameters will be performed elsewhere. 



Finite bias (slanting) is relevant to the situation of strong pinning of a dislocation at its ends in different 
vallyes of Peierls potential. Otherwise, dislocation will lower its energy by settling into a single Peierls 
valley. Strong end pinning with slanting imposes geometrical constraint and, therefore, introduces kinks 
even at T = 0. Such kinks change the responses dramatically: dislocation remains rough at T much lower 
than that in the absence of bias. Fig Q] demonstrates this effect in R2 response of non-slanted dislocation 
(lower blue curve) and of the slanted one (red upper curve). The first is smooth at low T and becomes 
thermally rough at elevated T. In contrast, the second, which has six geometrical kinks (given L x — 80, 
the bias is 6/80 = 0.075), exhibits just the opposite behavior - softening at low T. Accordingly, shear 
modulus G2, Eq. fTOl (see the inset in FigQJ hardens at low T for non-slanted dislocations and softens for 
slanted ones. Very recently such behavior has been reported in Refi^i for impurities free crystals of 4 Hc: 
shear modulus of 4 He crystals grown at 20mK from supcrfiuid has exhibited significantly smaller values 
than that at high T. Such softening, however, vanishes after the crystal is warmed up and cooled down 
again. According to RefJ^, this happens due to 3 He reentering solid 4 He and binding dislocations. 

We attribute such drastically different behavior of slanted dislocations in impurities free crystal to 



B. Finite Bias and shear modulus softening at low T 
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FIG. 4: (Color online) Spatial correlator Ri(x) (related to Ri as Ri — Y2 X Ri(x)) of the biased string with the 
same parameters as in FigQ] at different T . Spatial modulations reflect Wigner-type crystallization of the extra 
kinks. At high T, the spatial profile becomes that of the Granato-Liicke string pinned at both ends^. 



quantum nature of the extra kinks. Loosely speaking they form " superfluid" which enhances dislocation 
mobility at low T. The long-range potential ([5]) does not change significantly this picture in spite of 
inducing Wigner-type crystal quasi-ordcr among extra kinks. Such order is seen in the spatial oscillations 
of the correlation function R\{x) = (y(x)) / (RqF), FigQ] In Id achieving true crystal order is impossible, 
and, therefore, kinks remain "superfluid" and "solid" at the same time - some sort of "supersolid" of the 
extra kinks. 

However, dislocation network resides in 3d sample, and actual Wigncr crystallization of kinks should 
be possible. This implies that shear modulus may, first, demonstrate softening as T lowers and, then, 
at T < T* , - hardening and eventually the recovery to its ideal value G e . If, however, nd is very low, 
such recovery may never be observed. In order to estimate T* we compare it with typical "Coulomb" 
interaction cx Uc/r between kinks at their average separation r cx 1/y/Nd, provided their typical zero- 
point kinetic energy is smaller than potential. This gives T* ~ Ti^n^, where T\ stands for a typical 
energy scale 10K in 4 Hc. For dislocation densities Nd ~ 10 10 — 10 12 m~ 2 (the dimcnsionlcss density 
fid ~ i0~ 9 — ICP 7 ) we find T* ~ 1 — lOmK. This should be compared with the temperature of thermal 
roughening T r determined by typical energy of kink-antikink pair creation, which is a single dislocation 
effect and therefore is independent of rid- As discussed in Ref^, T r ~ 100 — 200mK. 



IV. SUPERCLIMB 



In the experiment^ pressure applied to superfluid 4 He inside the vycor "electrode" was transferred to 
the solid 4 He through relatively small contact area. Yet, this caused large uniform (which is inconsistent 
with plastic flow) accumulation of matter inside the solid. In this sense the effect of pressure can be viewed 
as a mean to change chemical potential \i. This phenomenon was called giant isochoric compressibility^ - 
macroscopic accumulation of 4 He atoms in response on changing [i. It has been proposed that the actual 
matter accumulation occurs due to building half-planes forming superfluid edge dislocations. 
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FIG. 5: (Color online) Temperature dependencies of K2 (with X2 oc R2, Eg. (|12[) 'l for superclimbing dislocation for 
different biases y m ax/L x ; F — 0.001, L x — 60, po(0) = 1, To = 0.2, K = 1, Uc = 1. Dot-dashed line represents 
renormalized superfluid stiffness p s {T), with the bare one taken in the form po(T) = po(0)(l — T/To) v , v = 0.67. 
Anomalously large response is found for one extra jog, bias=l/60 ~ 0.0167. 




FIG. 6: (Color online) Temperature and size dependencies of k(T,L x ) for superclimbing dislocation at zero bias. 
Size dependence of the peak (T = 0.05) can well be fit by the power law k(T = 0.05, L^) oc L™ , w m 1.32 
for Uc = 1. Inset shows p s (T,L x ). It becomes suppressed drastically at about 1/2 of To = 0.2 - mean field 
temperature where microscopic superfluidity vanishes. Size dependence of p s is insignificant if compared with k. 



We consider two isochoric compressibilities 
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AiV 



dN 



(11) 



Xi = 



NAfj,' 



Xi = 



Ndfj, 



where AA*" is the amount of He atoms accumulated in the sample upon applying chemical potential 
difference A[i, with N being total number of atoms in the sample. Variation of fi can be identified 
with deviation of local pressure p from the equilibrium value p as A/i = (p — po)Q where f2 stands 
for volume occupied by one atom. Force F in Eq. (J3j) is given by F = (p — po)b = Ajib/Vl. Thus, 
d.../dfi = (b/fl)d.../dF. AN can be related to average (dimensionless) density of climbing dislocations 
rid- Indeed, if the dislocation forest is characterized by typical length L x of smallest free segment, its 
displacement y changes number of atoms by ~ 3yLb per each block. Thus, the fraction of atoms changed 
is w 3ybL x /L x w ridy/b. This gives Xi,2 oc n^Ri^- For free segments of dislocations the responses 
become Rq oc L x which implies compressibilities comparable to liquid 4 Hc xo (since n c i oc 1/L X ). Thus, 
introducing relative compressibilities X1.2 — XL2/X0 an d R1.2 = Rx,2/Ro, one finds in the dimensionless 
units 



that is, xi,2 oc Ri^, where we ignore compressibility due to direct application of pressure to the area 
of contact between "electrodes' — and solid 4 He sample. The parameter 7 ~ 0.1 — 1 is the geometrical 
factor depending on the structure of the dislocation forest. 

The responses R\ t 2 have been determined from Monte Carlo simulations of the full action US]). In this 
paper we consider limit of small F only (where two responses coincide) . Non-linearities in the superclimb 
effect will be addressed elsewhere. As discussed in Ref4, Peierls potential with respect to jog formation 
is always relevant at low T, and, therefore, the compressibility will be suppressed. This is demonstrated 
in FigfS] where the response R2 becomes significantly suppressed at low T for, practically, any bias. As 
T rises, jogs are created thermally, provided the core superflow can still exist, making Peierls potential 
less and less relevant. Accordingly, Ri :2 should approach unity (that is, the free string limit i?o oc L x ). 
At higher T superfluidity ceases to exist and so does the giant compressibility effect. 

As discussed in RefJ^, dislocation network superfluidity is characterized by two temperature scales - 
microscopic To ~ IK (comparable with the bulk A-tempcrature) and macroscopic T c oc Tqci/L x << Tq. 
In our model, we choose the bare superfluid stiffness in Eq.® to vanish according to the 3D XY- model 
law po(T) = po(0)(l — T/Tq) u , v « 0.67 with the choice Tq — 0.2 (in dimensionless units). As can be 
seen in Fig [51 high-T behavior of R2 reflects suppression of the renormalized superfluid stiffness p s (T) 
which happens at about 1/2 of To. Strictly speaking, the 3d form of the bare stiffness is not valid very 
close to To and one should rather use the mean field form. However, since major suppression occurs at 
much lower T, the actual form of po(T) is not important. 

Another prediction for the superclimb effect is that the core superfluidity becomes of non-Luttinger 
type with the free-particle- type dispersion^ for Uc = 0. In terms of the superfluid compressibility k this 
implies strong dependence k oc L x on the dislocation length in the temperature range where superclimb 
is observed. We have measured superfluid compressibility through imaginary time windings fluctuations 
of the J-currents m t . Fig(5] presents temperature dependence of n for dislocation lengths L x = 10 — 100 
at zero bias. At low T all curves collapse, practically, into a single one. This indicates the region of 
the standard Luttinger liquid behavior. As T rises, curves for different L x diverge, with their maxima 
obeying asymptotic behavior n max oc L™, w ~ 1.32 for Uc = 1. In general, the exponent w depends 
on long-range interaction strength Uc, approaching w — 2 for Uc = 0. At higher T the curves collapse 
again and eventually k vanishes. 

We have found that bias with respect to climb does not render Peierls potential irrelevant at T = 0. 
It simply lowers the temperature when xi,2 become suppressed. As can be seen from Figj5j \2 a t finite 
biases starts diminishing with lowering T at approximately 3 times smaller temperature than that for 
zero bias. Qualitatively, this scenario has been found to be independent of the "Coulomb" interaction 
([5J between jogs. 



Xl,2 = 7^1,2, 



(12) 



V. CONCLUSIONS 
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Low-T dislocation glide is characterized by hysteretic behavior and ratchet-type shape of i?i-response as 
a function of applied stress, with the i?2-response exhibiting peaks at the steepest parts of the "ratchet" . 
The threshold stress is determined by Peierls gap and it vanishes inversely proportional to approximately 
dislocation length. Thus, these type of measurements may reveal temperature dependence of Peierls gap 
as well as give an estimate of a typical length of dislocation segments between strongly pinned crosslinking 
points^. 

Geometrical bias (due to ends of dislocation residing in different Peierls potential valleys) introduces 
kinks even at T = and produces dramatic effect on gliding dislocation: while at zero bias dislocation 
is always quantum smooth and shows thermal roughening above some temperature Tji 3 -, single biased 
dislocation remains rough at T << T r . Furthcmorc, effective rigidity of such dislocation shows softening 
with respect to its value at high T. If such biased dislocations dominate the network, shear modulus 
should exhibit softening at low T— - the effect opposite to the observation^. Long-range forces between 
kinks lead to Wigner-type quasi-crystalization of kinks. This effect, however, does not affect significantly 
quantum roughening of biased dislocation due to its Id nature. At finite density of dislocations, true 
Wigncr crystalization of kinks should take place below temperatures determined by dislocation density 
when the system becomes essentially 3d. Such softening can be suppressed by 3 Hc impurities localizing 
the extra kinks. Thus, it is important to study shear moduli in impurities free solid 4 He crystals as 
suggested in Ref.— in order to reveal the nature of dislocation crosslinking. Observing low-T moduli 
softening insensitive to thermal cycling would imply significant presence of slanted dislocations stabilized 
by strong pinning at the crosslinking. Conversely, observing hardening induced by annealing would mean 
that the dislocation forest is able to reach more stable configuration without geometrically imposed zero- 
point kinks. As recent report^ indicates, despite restoring the low-T hardening after thermal cycling, 
there is very slow relaxation toward the low-T soft state due to escape of 3 He back into the liquid which 
is in a contact with the solid. This favors the first scenario - stable dislocation network dominated by 
slanted dislocations. 

Isochoric compressibility due to superclimbing dislocations has been evaluated and the prediction of 
its low-T suppression^ has been confirmed. In contrast to glide, such suppression is not qualitatively 
affected by geometrical slanting. We attribute this to the nature of jog superclimbing dynamics - strong 
control of jog displacement by supercurrents. More specifically, in order to displace Nj jogs by some 
distance S along the dislocation, SNj atoms must be transferred through its ends. Two factors inhibit 
such transfer - interaction between (extra) jogs and low-T suppression of compressibility due to Peierls 
barrier far from jogs (where bias is irrelevant). Single jog is not affected by the interaction and, thus, can 
move freely to induce large fluctuations of the total number of atoms, that is, - large compressibility. At 
low T the second factor suppresses this anomalous single-jog behavior. An interesting subject for future 
studies is non-linearities of the superclimb with respect to applied pressure. 
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